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We propose an easy-to-build easy-to-detect scheme for realizing Majorana fermions at the ends 
of a chain of magnetic atoms on the surface of a superconductor. Model calculations show that 
such chains can be easily tuned between trivial and topological ground state. In the latter, spatial 
resolved spectroscopy can be used to probe the Majorana fermion end states. Decoupled Majorana 
bound states can form even in short magnetic chains consisting of only tens of atoms. We propose 
scanning tunneling microscopy as the ideal technique to fabricate such systems and probe their 
topological properties. 
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The interest in topological quantum computing and 
non-abelian braiding has inspired many recent propos- 
als to create Majorana fermions (MFs) in various experi- 
mental systems. Following Kitaev's seminal proposal [1], 
many approaches have been considered including those 
based on topological insulators [2, 3]; atoms trapped in 
optical lattices [4-6]; semiconductors with strong spin- 
orbit interaction in two and one dimension [7-9]; cou- 
pled quantum dots [10, 11]; and those that combine mag- 
netism of and superconductivity [12-14]. The aim of 
these approaches is to create a topological superconduc- 
tor in which MFs emerge as the single excitations at the 
boundaries. Since MFs are its own antiparticles, they 
are predicted to appear in tunneling spectroscopy exper- 
iments as zero bias peaks [15-18]. Such peaks have been 
indeed observed in several experiments and interpreted 
as the signatures of MFs [19-21]. However, these experi- 
ments are not spatially resolved to detect the position of 
the MFs. Additionally, in many instances, the presence of 
disorder can result in spurious zero bias anomalies even 
when the system is not topological [22, 23]. It is therefore 
desirable to identify easy-to-fabricate condensed matter 
systems in which MF can be spatially resolved and dis- 
tinguishable from spurious disorder effects. 

In this letter, we theoretically investigate conditions 
for which a chain of magnetic atoms on the surface of an 
s-wave superconductor can host MF modes. We explore 
the parameter space for which this system is topological 
and show that even relatively short chains made of only 
~ 50 atoms can host robust localized MFs. Our proposed 
structures can be fabricated using scanning tunneling mi- 
croscopy (STM), which has previously been used to as- 
semble structures of various shapes with tens of atoms us- 
ing lateral atomic manipulation techniques [24-26]. Spa- 
tially resolved STM spectroscopy of such disorder-free 
chains can be used to probe the presence of MF end 
modes. 

As shown in Fig. 1, we consider an array of magnetic 
atoms (such as atoms of 3d or 4f metals with a net mag- 
netic moment) which are deposited on a single crystal 



surface of an s-wave superconductor (such as niobium 
(Nb) or lead (Pb)) and arranged into chains using the 
STM. The interaction of a single magnetic moment with 
the superconductor gives rise to the so-called Yu-Shiba- 
Rusinov states [27-30] that have been previously detected 
from both 3d and 4f atoms on the surface of Nb and Pb 
using an STM [31, 32]. The results of these previous ex- 
periments (with Gd and Mn deposited on Nb) agree well 
with model calculations in which the magnetic moment 
is assumed to be static [31, 33, 34]. In addition, recent 
spin polarized STM studies indicate that in magnetic ar- 
rays with > 10 atoms spin dynamics is greatly suppressed 
[35] . It is therefore reasonable to model moments of mag- 
netic atoms as static classical spins. In general, magnetic 
moments in these chains can form various configurations 
including a spiral [36] . 




Figure 1. Schematic of the experimental setup. An array 
of magnetic atoms (red spheres) is assembled using scanning 
tunneling microscope on the surface of s-wave superconduc- 
tor (gray background). The system is modeled by the two 
dimensional N a x Nb array in which magnetic atoms are em- 
bedded (inset). Throughout the paper we consider the case 
where magnetic moments are in the plane defined by N a and 
Z direction. 



To describe this system we use a two-dimensional tight- 
binding model Hamiltonian of an s-wave superconductor 
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with an array of magnetic atoms : 
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The operators / and /' correspond to electron annihila- 
tion and creation respectively, t is the hopping amplitude 
between adjacent sites <i,j> of a two-dimensional lattice, 
fj, is the chemical potential, Aj is the local superconduct- 
ing gap associated with a host superconductor (equal to 
Ao in the absence of magnetic atoms) . The effective mag- 
netic field B n gives rise to a local Zeeman energy on the 
atoms which are arranged in a one- dimensional array of 
sites {n}. We consider the case of identical atoms, i.e. 
\B n \ = B. Throughout the paper we normalize all simu- 
lation parameters to the value of A . 

In order to obtain the two-dimensional gap profile in 
the vicinity of the atomic chain, we self-consistently solve 
the resulting Bogoliubov-de Gennes equations (BdG) 
[37]. We assume a constant on-site pairing coupling V 
for a grid of iV a x Nf, lattice sites in the middle of which 
N a local magnetic moments with strength B are embed- 
ded (see [38] section 1 for details). The calculations are 
performed with open boundary conditions (BC) in the 
Ni, direction, and both open and periodic BC in the N a 
direction to show the presence or absence of MF at the 
end of the chain and to compute the Pfaffian index (Pf ) 
[1[. Previous calculations showed that a single magnetic 
moment gives rise to a state inside the superconducting 
gap that has an energy close to Ao for low B. As the 
value of B is increased the energy of this state is contin- 
uously tuned to zero [33, 34, 39]. This zero crossing is a 
signature of a quantum phase transition, at which the im- 
purity site traps a single quasi-particle [28, 40]. A similar 
phase transition occurs in the case of a few magnetic mo- 
ments [39, 41]. The transition obviously coincides with 
a change of the sign of the Pfaffian (computed in a peri- 
odic geometry) for the system, indicating a change of the 
fermion parity in the ground state. This is the charac- 
teristic signature of a topological non-trivial phase with 
MF end modes [1[. 

An example of a transition into a topologically non- 
trivial phase for our atomic chain is illustrated in Fig. 
2, which shows the lowest energy level and the Pfaffian 
as a function of B in the case of 96 magnetic moments. 
The angle between adjacent magnetic moments, 0, plays 
a key role in determining whether this system is topo- 
logical (see below), and has been assumed to be 27r/3 
for the results shown in Fig. 2. The most important 
feature of this calculation is that in the parameter win- 
dow 2.2 < B/Ao < 3.45, in which for periodic BC in 
the iV a direction the Pfaffian is negative, and the spa- 
tial extent of the lowest excited state (Fig. 2b) (for open 
BC) shows the presence of MFs at the ends of the chain. 



This behavior can be contrasted with that of B/Aq = 2.1 
(Fig. 2c). In this case Pfaffian is positive and the low- 
est energy excitation is distributed approximately evenly 
along the chain. A calculation of the local density of 
states (LDOS) as a function of energy shown in Fig. 2d 
clearly demonstrates that the topological case shows a 
zero bias peak associated with MF when tunneling at the 
end of the chain, while the middle of the system exhibits 
a mini-gap. In the non-topological phase sufficiently far 
away from the transition point, the system shows a clear 
gap throughout the chain and absence of zero energy end 
modes (Fig. 2e). 

The emerging MF end modes considered here are lo- 
calized on a very short length scale at the last few sites 
of the atomic chain. This situation can be contrasted to 
the proposals involving semiconductor nanowires in prox- 
imity with superconductors, where the coherence length 
of the superconductor sets the length scale for MFs [9[. 
The spatial extent of our MFs is reminiscent of the extent 
of the Yu-Shiba-Rusinov states created by single atoms, 
which have been shown both experimentally and theoret- 
ically to decay on length scales associated with the Fermi 
wavelength of a superconductor [31, 34]. Note that these 
states do have long tails associated with the supercon- 
ducting coherence length, however this decay is strongly 
enhanced with an algebraic decay pre-factor [33, 34]. 

While we used a self-consistent BdG calculation for 
realistic modeling of experimental situation, a more effi- 
cient approach to gain physical insight into this system 
is to consider an effective ID model of magnetic atoms 
on superconducting sites, which is just the Nf, — 1 limit 
of our 2D model. Note that in ID, all information about 
the superconductor is simply included in the strength 
of the on-site s-wave gap A and the hopping term de- 
scribes coupling between the impurities on superconduct- 
ing sites only (as opposed to the superconductor band- 
width in BdG model above, see [38] section 2). Fig. 3 
shows that a ID model qualitatively gives similar results 
the 2D model. Importantly, the hopping term, which 
can be tuned experimentally by placing atoms at differ- 
ent distances, may also drive quantum phase transition 
from the trivial phase (Pf>0) to the topological phase 
(Pf<0) with MFs at the ends. A one-dimensional ver- 
sion of this Hamiltonian is also considered in Ref. [12] in 
the context of MFs in disordered magnetic islands on a 
superconductor . 

A key advantage of the ID model is that it lends itself 
to an analytical solution, which shows that for a given 
angle 9 between adjacent moments, the Pfaffian for the 
system is negative when 

v /A2 + (H+2|tcos(0/2)|)2 > \B\, 

\B\ > ^At + (\,4-2\tcos(8/2)\y (2) 

(see [38] , section 3 for the derivation) . The negative value 
of the Pfaffian is a necessary condition for this system to 



3 




Figure 2. (a) Calculated energy spectrum, marked by 
blue lines, for 96 classical spins placed in the middle of the 
N a x Nb — 96 x 19 grid using periodic BC. Parameters for 
the plot are: /i/A = 2.12, t/A = 2.34, V/A = 2.81 and 
T/Ao = 0.01. The regions corresponding to the trivial phase 
(Pf>0) are shaded gray. Red thick line represents the lowest 
energy excitation using open BC. (b,c) The spatial distribu- 
tion of the local density of states corresponding to the lowest 
excitation state in the non-trivial (B/Ao = 2.87, Pf < 0) and 
the trivial (B/ Ao= 2.23, Pf >0) phase. Lattice coordinates 
X and Y correspond to the N a direction (along the chain) 
and Nb direction (orthogonal to the chain) respectively. (d,e) 
Local density of states at the chain ends (blue solid line) and 
in the middle of the chain (gray dashed line) as a function of 
energy for non-trivial and trivial phase taking into account 
first 96 energy eigenvalues. The intrinsic line- width of the 
energy eigenstates is taken to be u/Ao = 1 x 10 -3 for this 
plot. 



be in a topological phase; however, it not sufficient, as 
the bulk of atomic chain remains must also be gapped. 
For example, 8 = 0,tt have the widest range of negative 
Pfaffian in Eq. 2; unfortunately, this full range is gap- 
less. The min-gap for low energy excitation is related to 
strength of the p-wave pairing that emerges on the chain 
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Figure 3. (a,b) The spatial profile of the two lowest excitation 
states of magnetic chain containing 48 atoms for /x/Ao = 4, 
B/Ao = 5, 6 = 7r/2 (marked by red solid and green dashed 
line respectively) . Tuning the hopping term t drives quantum 
phase transition from the trivial (i/Ao = 0.4, Pf > 0) (a) to 
the topological (t/Ao = 1, Pf <0) phase (b). (c,d), Local 
density of states calculated for the same parameters as in (a) 
and (b) at the chain ends (blue solid line) and in the middle 
of the chain (gray dashed line). Note that for this choice 
of parameters spectrum in (c) is asymmetric in energy (see 
inset). Importantly, in (d) the two MF states around zero 
energy are separated by the effective mini gap A p from the 
other states in the spectrum (marked by double arrow line). 



because of the combination of hopping, pairing, and lo- 
cal Zeeman terms in the Hamiltonian. Calculations of 
the spectra in both 2D and ID model described above 
reveal the energy scale, which separates the zero energy 
MF states (localized at the two ends) from the next avail- 
able excitation of the system. In a certain limit, the ID 
model can be directly mapped [12] to the original pro- 
posal by Kitaev for realization of MF end mode, which 
is a superconducting wire with nearest neighbor pairing 
[1], but general eigenvalues can be obtained even without 
this mapping, see [38] section 2. The value of this mini- 
gap depends on the relative values of fi. t, B, and angle 
9 (see Fig 4) . 

A non-collinear arrangement of magnetic moments in a 
chain is essential to realize robust MF end modes. When 
transformed to a basis parallel to the spiraling on-site 
magnetic field, the hopping becomes spin-dependent giv- 
ing rise to spin-orbit coupling and hence to the usual 
mechanisms for MF end modes. Without detailed mod- 
eling of the surface magnetism it is difficult to predict 
whether specific magnetic atomic chains would have a 
spiral spin-arrangement. We suggest that exploring the 
full freedom of the linear chain geometry may provide 
a feasible approach to create favorable conditions for 
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Figure 4. The value of the mini gap as a function of tunnel 
coupling and angle 8 calculated for the ID model : (a) fi/Ao = 
2 and B/A = 3; (b) n/A = 2 and B/A = 5 ; (c) fi/A = 5 
and B/A = 2; (d) fi/A = 5 and B/A = 4. 

non-collinear magnetic moments of adjacent atoms. For 
example, double or zig-zag chain structures with anti- 
ferromagnetic interactions are likely to become frustrated 
and result in spiral orientation of magnetic moments in 
the chain [36] . To explore some of these possible geome- 
tries (Fig. 5a), we map these chains into equivalent linear 
chains with the nearest t\ and the next nearest i 2 hop- 
ping as shown in Fig. 5b. In the simplest case for which 
9 is assumed constant, we show that these chains can also 
support topological phase when 

VA 2 + (/i + 2 cos(0/2)ii - 2 cos(6%) 2 > \B\, 
\B\ > yfA 2 + (/i - 2 cos(0/2)£i - 2 cos(6%) 2 (3) 

(see [38] section 3a for further details). We note again 
that these chains may provide easy-to-fabricate struc- 
tures that would ensure non-collinear spin arrangements 
required for realization of MF end modes. 



ductor (Fig. 4). Nevertheless, using an s-wave supercon- 
ductor with large gap Ao (and measuring at the lowest 
temperatures) would increase the chance of experimen- 
tal success. Other factors such as size of the magnetic 
moment B or hopping matrix element t are also impor- 
tant and can be optimized experimentally using magnetic 
atoms with different spin or building chains with different 
spacing. A systematic experimental approach can start 
by characterizing the single-impurity states and their 
modification when impurities are brought close enough to 
interact [32] . These measurements could be used to map 
effective ID model parameters (effective hopping, chemi- 
cal potential and exchange coupling) and allow investiga- 
tion of the finite size effects on the excitation spectrum. 
A different approach would be to start from magnetic 
chains grown using self-assembled techniques. Note that 
self-assembled chains consisting of ~ 50 atoms with spiral 
arrangement of magnetic moments are already reported 
[36]. Such chains would be an ideal starting point to 
investigate interaction between Majorana fermions. For 
example, examining coupled chains can provide direct ex- 
perimental means to demonstrate the Z2 character of the 
MF end modes by showing that they appear only in odd 
number of coupled chains. Finally, as structures of differ- 
ent shapes are equally easy to assemble in STM, one can 
envision viable route towards braiding experiments in ar- 
rays of coupled chains in a similar fashion as proposed for 
semiconductor nanowire structures [42-44]. 
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Figure 5. (a) Array of magnetic atoms arranged in two rows 
(zig-zag chain). The coupling among neighboring atoms cor- 
responding to different rows is t\ and the coupling between 
atoms within the same row is t^. (b) Equivalent magnetic 
moment configuration represented as a single chain with the 
next nearest coupling. 

Lastly, we comment on the experimental feasibility 
of the proposed approach. As shown here the strength 
of the mini-gap associated with the p-wave pairing can 
sometimes exceed 30-40% of the gap of the host supercon- 
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1. DETAILS OF THE TWO-DIMENSIONAL MODEL 

The numerical calculations of the two dimensional model are preformed self-consistently [SI]. First we solve Bogoli- 
ubov de Gennes equation for the Hamiltonian of our system (Eq. 1) and find electron- and hole- like wave functions 
u na(i,j) and v ncr (i,j). The excitation spectrum is then used to determine local on-site value of superconducting gap 

A(i,j) = -V/2^T{(2/ n - l)(u n ,t(i,j)X;(i.i) + «n4(i»jX,t(i»J'))} ( S1 ) 

n 

where V is the pairing strength, /„ = l/(e e "/ T + 1) Fermi function, n sums over energy eigenvalues e„ and T is the 
temperature. In the next iteration the obtained values A(i,j) are used to recalculate the excitations spectrum. This 
iterative process is stopped after 20 iterations, typically assuring relative error of A(i,j) to be less than 10 -2 . In the 
last step we also calculate the spatial profile of local density of states (LDOS) defined as 

p(h.j, e ) = ^2(\ u n,a(i,j)\ 2 S{e - e„) + \v nttr (i, j)\ 2 S(e + e n ). 

n,a 

We used open boundary conditions for the Nb direction and periodic and open boundary conditions in the N a direction. 
Periodic boundary conditions are used to compute Pfaffian index of the system (as discussed in the section 2). When 
the Pfaffian is negative (i.e. in topological phase) open boundary conditions are used to verify the existence of the 
two zero energy modes with LDOS localized at the chain ends. These two states are identified as a pair of Majorana 
fermions (MFs). 

Although self-consistent calculations are useful to explicitly verify the occurrence of MFs in more realistic conditions, 
they require intensive computational resources. Qualitatively similar results can be obtained in a toy model without 
imposing self-consistency. This type of calculations allow to quickly investigate the parameter space for which the 
system is in topological phase. Further, in the one-dimensional case (for AT& = 1) analytical conditions can be found 
in some cases as shown in the following sections. 



2. HAMILTONIAN IN ONE DIMENSION AND THE TOPOLOGICAL INDEX 

In the case of N b = 1, the Hamiltonian in Eq. 1 reduces to: 

H = ^ tnfl a fn+la + t*nfl+lafna ~ fi ^ f la fna + ^{^n ' v)a(lflafnP + ^2 Ao -^T^ + A o/"4-/nt ( S2 ) 
na na na/3 n 

In this one-dimensional Hamiltonian we assume a fixed value of Ao for every site. As noted in the main text, for this 
model, tunnel coupling t defines solely hopping term between magnetic atoms. Note that this is an experimentally 
tunable parameter since it can be changed by placing atoms on different distances. To define Pfaffian of this system, 
we first rewrite the Hamiltonian using the real-fermion Majorana representation. The mapping from complex fermions 
to MF basis can be done in the following way: 

Q>2n— 1,(T fna ~t~ fna' &2n,a — ^{fna fna)- a tia — &na i {^nffi^ma'} 2S nm S a a i . (S3) 

Note that the site index of the complex fcrmion / goes from 1 to N a while the index of the real Majorana fermions 
goes from 1 to 2N a . Since for each complex fcrmion f na there is two real MFs a,2na and a 2n -ia the number of sites 
is effectively doubled. In Majorana representation ID Hamiltonian reads as 

tj v~ *N iRe(t n ) j \ . ilm(t n ) / . \ 

tl — Z^n=l 2 ( a 2n-laCl2n+2a ~ ^na^n+la) H g l a 2ri-l<ra2n+l(T + 0,2nad2n+2a ) — 

— l -2 0-2n-laO-2na + (a2ri-lt a 2rit — a 2n-llO<2nl) + l -^f 3L (a2ri-lt a 2ni _ «2nt a 2n-li) 

_ i^tS.(a 2 n-lt a 2n-14. + a 2rit«2r4) + ^ {d-2n- H a 2nt + a 2r4a2ri-lt) (S4) 
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We have assumed a real order parameter A = Aq. Importantly, Hamiltonian in this representation can be rewritten 

as 

i 2N 

H = - a l<jMa:ma'a m <y' (S5) 

/,m— 1 

where A\ a . m(J , = A la ■ ma r — — A m( yi.i a is an AN a x AN a (2N a MFs per spin species) antisymmetric matrix with elements 
easily readable from the above expression. For the upper defined form, we define the topological index as 

sign(P. fC-A la[ma ,)) (S6) 

which effectively describes the topological phase of the system [S2]. Any change between +1 and —1 of the index 
indicates a phase transition where the system changes its topological phase. Note that the extension of this index 
to the 2D model is straightforward and was done using the same majorana representation. Starting from a (2 x 
N a x Nb) x (2 x N a x iVft) matrix representation of the Hamiltonian (Eq. 1) in the complex fermion basis we get 
(4 x N a x Nb) x (4 x N a x Nb) Hamiltonan matrix in Majorana basis. The corresponding Ai a . m<7 ' matrix is again 
antisymmetric and the same definition of the topological index can be used. 



3. ANALYTICAL SOLUTION FOR TOPOLOGICAL PHASE IN THE CASE OF SPIN HELIX 

In this section we solve the one-dimensional Hamiltonian (Eq. S2) for periodic boundary conditions when angle 
between adjacent magnetic atoms is constant. Since all atoms are identical the magnetic moments on each site have 
same magnitude but differ in direction: B n = B n (sm(9 n ) cos(0„)x + sin(0 n ) sm((/> n )y + cos(6 n )z). We start from the 
complex fermion representation and first rotate the fermionic operators into an t-i basis on each atomic site using 
unitary rotation described in Rcf. [S3] 

fnt \_rrf 9n t \_f COs(0„/2) - S hl(0„/2) e -^ \ / 9nt \ 

U J' n \ 9ni ) \ sin(0„/2)e^ cos(0„/2) ) \ 9ni ) 

As U n is a unitary transformation fermions g na have the same anticommutation properties as f na . This transformation 
doesn't change the form of the chemical potential and gap terms, whereas the magnetic moment term becomes 
diagonal. The hopping term now acquires inter-spin components, and is no longer diagonal in spin. The rotated 
Hamiltonian now reads as 

H = tnttn^pgladn+lfi + *^,/},a0l+lc«5n/3 + B (T za pgl a g nfj - fl ^ g\ a9 n a + A (.9^5^ + 9ni9n^) (S8) 

n,a,f) na n 

where 0„ is defined as 

0„ = UtU n+1 = ( 2 -£ ) (S9) 

with matrix elements being a n — cos(0„/2) cos(0„ + i/2) + sin(#„/2) sin(#„ + i/2)e _ ^™~^™+ 1 ) and j3 n = 
— sin(#„/2) cos(0 n +i/2)e l< ^™ + cos(6>„/2) sin(0„ + i/2)e^™ +1 . In the case of magnetic moments aligned in the x-z plane 
with constant angle 9 between adjacent moments coefficients a n and /?„ reduce to a = cos(#/2) and j3 — s'm(9/2) 
respectively. We rewrite the rotated operators g na in Majorana basis as follows: 

9n<7 = ^(&2n-l<r + ifona), J?L = ^(&2n-la - ifona) ■ (S10) 

The b Majorana operators are the rotated a Majorana operators in the basis parallel to the magnetic field. In this 
basis, the Hamiltonian becomes: 

H = tn ^ an i (&2Tt-lt^2n+lt + &2rtt&2n+2f _ ^n-H^n+H _ &2r4&2n+2.|.) + 

+ ''"^f""^ (^2n-lt^2w+2t — b 2n -fb 2n +l-f + b 2n -llb 2 n+2l — b 2n \b 2n+ il) + 
+ En ^\ ^ {b 2 n- l-fb 2 n+ll + b 2n -\b 2n+2 ^ _ &2n-l.|.&2n+lT — ^2r4^2n+2t) + 

_l it n (Pn+P n ) (b 2n _ 1 ^b2n+2i — b 2n -fb 2n+ i_i + b 2n -i^b 2n+2 ^ + b 2n ^b 2n+ i-f) + 

+ En ^(fon-lt^nt - kin-l^ni) - ^ (&2n- lf^nf + fan-l&nl) + ^ (&2n-l.L&2nt - &2n-lf&2n4.) (Sll) 
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To explicitly calculate the Pfaffian index and energy eigenvalues we need to Fourier transform MF operators as 

fcn-i* = 4= E e~ iqn K^, b 2na = -L e~ iqn b q ,2,* (S12) 

v n v n 

where 1, 2 are the two flavors of MF operators coming from splitting the on-site fermion, and q is the Fourrier vector 
which takes values to be determined below (there are several cases for the values of q depending on the boundary 
conditions, which in turn depend on the value of the rotation of the spin). We can then write the Hamiltonian in the 
basis b q = (b g ,i,t, fyz,i4> b q,2^, 

H =~Jl b l A ^)b q (S13) 

Q 

where A(q) matrix has nonzero matrix elements: 

A 13 (q) = 2ta cos(q) + B - fi, A 14 (q) = 2it(3 sin(g) - A (S14) 



-423(g) = -2it/3sm(q) + A A 24 (q) = 2tacos(q) - B - fj, (S15) 

Periodic boundary conditions restrain the angle 9 to be 9 — jf-Nh where Nh defines number of times the magnetic 
moment is rotating in an N a site system. In other words, the magnetic moment rotates by an angle 2nNh going from 
site n = 1 to site n — N a . A further complication may arise since we are rotating a spin 1 /2 system. If the magnetic 
moment rotates an odd number of times, in the rotated basis we have a situation where spin 1/2 is rotating an odd 
number of times giving rise to a —1 sign. This changes the periodic boundary conditions to anti-periodic. Therefore, 
the allowed values of the Fourier vector q depend on whether N h and N a are odd or even. The four cases with allowed 
values of q are summarized as follows: 

If mod [N h , N a ] < N a /2 and if N h even 

If mod [N h ,N a ] <N a /2 and if N h odd 

If mod [N h ,N a ] > N a /2 and if mod [N a - mod [JV h ,iV ],2] =0 [ ' 

If mod[N h ,N a ]> N a /2 and if mod [N a - mod [N h , N a ], 2] = 1 

Replacing the allowed q in the expression for Pfaffian one can obtain topological index: 

Sign(P f[A(q - 0)])Sign(P/[A(g = tt)]) = Sign([B 2 - (A% + (/i + 2atf)}[B 2 - (A„ + 2atf)}). (S17) 

Note that q = 0, tt values are allowed only in only when Nh and N s are both even. The nontrivial topological phase 
appears when the above index is —1 i.e. when: 

^/A§ + (H+2H)2 > |B | > ^/A§ + (H-2H)2. (S18) 

This is the necessary condition for Majorana modes to appear in the one dimensional model assuming constant 9. For 
actually having Majorana modes the system also has to be gapped. Therefore it is necessary to calculate the energy 
spectrum which can be easily obtained by diagonalizing the Eq. S13. For each q four eigenvalues are 

E{q) = ±^B 2 + 4(1 - a 2 )t 2 sm 2 (q) + {fi - 2atcos{q)f + A 2 ± 2<J B%A 2 + (p - 2at cos{q)) 2 (B% + 4(1 - a 2 )t 2 sin 2 (q)) 

(S19) 

where all possible combinations of the ± signs are possible. At q = 0, the bands closest to the energy become 
±(B - ^/Al + (fi- 2at) 2 ) and at q = n they are ±(B - ^ A^ + (fi + 2at) 2 ) . This shows that when the Pfaffian 
changes sign at the edge of the topological phase the gap also vanishes. 

Notice that for anti-ferromagnetic coupling (9 = n) the criteria in Eq. S18 collapses to a single point i.e for this 
magnetic moment orientation system never crosses into topological phase. Naively, one could expect that topological 
phase is maximized for ferromagnetic coupling (9 = 0). In this case the energy eigenvalues (Eq. S19) reduce to 
E(q) = ±B ± ^/Aq + (fi — 2tcos(q)) 2 . Unfortunately, for this case, in the whole range where Pfaffian is negative the 
system is also gapless and therefore it does not have topological phase. 

The absence of MF modes in the ferromagnetic and anti-ferromagnetic configurations can be also confirmed using 
a self consistent two dimensional model as illustrated by example shown in Fig. SI. In the anti- ferromagnetic 
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FIG. SI. Energy spectrum showing lowest 24 states of the magnetic chain for anti-ferromagnetic (a) and ferromagnetic (b) 
orientation of magnetic moments (N a x JVj, = 24 x 25, t/Ao = 3.83, /u/Ao = 3.19, T/Aq = 0.012). For anti-ferromagnetic 
coupling (a) the topological index stays positive (Pf > 0) in the whole range of B. For ferromagnetic coupling PfafBan changes 
sign around B/Aq m 6 but the mini-gap in this range (for B/Aq > 6) practically does not exist. Note that at around B/Ao ~ 5 
the lowest energy levels also come very close to zero energy, but the Pfafnan does not change sign around this point. 

configuration the levels approach zero energy but the two lowest levels never cross, (Fig. Sl(a)). Therefore, the chain 
remains in the trivial phase (Pf > 0) for the whole range of B. On the other hand, for the ferromagnetic configuration 
several levels practically collapse as the lowest energy level crosses zero, Fig. Sl(b). In this configuration Pf < but 
there is no well defined gap of the system. In other word system is gapless and the MF modes can not be therefore 
separated from the remaining states in the spectrum. 

3a. Topological phase of the simple zig-zag chain 

The analytical results above can be extended to more complex chains. In this subsection we consider zig-zag chains 
of impurities where atoms are placed in two rows as shown in Fig. 5 of the main text. This configuration is interesting 
since, for antiferrromagnetic coupling, it represents a frustrated spin system with potentially ground state where spins 
are canted off from purely anti-ferromagnetic alignment. A zig-zag chain can be thought of as a skewed ladder, where 
one of the chains in the ladder is shifted mid-bond with respect to the other chain. In this case, the zig-zag chain can 
be thought of as a single chain with both nearest neighbour hopping t\ n and second neighbour hopping t% n . To find the 
Hamiltonian for we can again project to the on-site spin basis B n = B n (sm(d n ) cos(<fi n )x + sm(9 n ) sm(<fi n )y + cos(9 n ) z) 
with B n > and the Hamiltonian becomes: 

H = *l™^l™,a,/3sLSn+l/3+ < 2n^2ri,a,05L.%+2/3+^ (S20) 

n,a,P na n 

where the matrix f2 lrl is identical to the one before while: 

n 2n = ulu n+2 = ( a2n ) (S2i) 

\ P2n Ct 2n J 

where a 2n = cos(0„/2) cos(#„ +2 /2) + sin(0 n /2) sm(9 n+2 /2)e~^- 6 ^ and f3 2n = - sin(0„/2) cos(0 n+2 /2)e i *» + 
cos(0 n /2) sm(d n+2 /2)e t9 "+\ 

If we now assume that the magnetic moment form spiral configuration, just like in the single-chain case, the situation 
is easy to analyze. We take t\ n — t\, t 2n — t 2 constant across the chain. By taking (j> n = and 8 n = (n — l)j^-Nh in 
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which the spin is rotating N h times in an N s site zig-zag chain, we have: 

a ln = ax = cos(n N h /N s ), ln = fa = sm(nN h /N s ), a 2n = a 2 = cos(2n N h /N s ), f3 2n = f3 2 = sm(2wN h /N s ) 

(S22) 

By going to the Majorana basis and Fourier space the Hamiltonian can again in the same form as Eq. S13 H = 
\ ^2 q b\A(q)b q with nonzero elements of the matrix A(q) being 

A 13 (q) = 2t iai cos(q) + 2t 2 a 2 cos(2<?) + B Q - fj,, A 14 (q) = 2it x fa sin(g) + 2it 2 /3 2 sin(2q) - A , (S23) 

A 2 z{q) = -2it t fa sin(g) - 2it 2 (3 2 sin(2<?) + A , A 24 (q) = 2t x a x cos(g) + 2t 2 a 2 cos(2g) - B - fi. (S24) 

This result is simple extension the introduction of the second neighbor gives just another added hopping contribution 
in the Majorana basis. The energy spectrum of the problem reads: 

E(q) = ±\^B~l + a{q) + b(q) +A%± 2^/B 2 A 2 + b(q)(B 2 + a(q)) (S25) 

where a (q) = (2fatism(q) + 2f3 2 t 2 sin(2q)) 2 , 6(g) = (p, — 2a.\t\ cos(q) — 2a 2 t 2 cos(2g)) 2 corresponds to replacing the 
2f3ts'm(q) — > 2fati sin(g) + 2(3 2 t 2 sin(2g) and 2at cos(g) — > 2t\ot\ cos(g) + 2i2C*2 cos{2q) in the simple chain problem. 

The topological index Pfafnan of the problem also comes out as a simple substitution but since the harmonics of the 
next nearest neighbor term is cos(2g), this term does not change sign between q = and q = n. The final condition 
for nontrivial phase is given by following expression 

^ 'a 2 + (ju + 2aiti - 2a 2 t 2 ) 2 > \B \ > A 2 + (fi - 2 ai h -2a 2 t 2 ) 2 . (S26) 

Note that the existence of a next nearest coupling, in addition to possibly providing favorable magnetic configuration, 
effectively shifts the chemical potential by — 20:2^2 • 



[SI] P. D. Sacramento, V. K. Dugaev, and V. R. Vieira, Phys. Rev. B 76, 014512 (Jul 2007) 
[S2] A. Y. Kitaev, Physics-Uspekhi 44, 131 (2001) 

[S3] T.-P. Choy, J. M. Edge, A. R. Akhmerov, and C. W. J. Beenakker, Phys. Rev. B 84, 195442 (Nov 2011) 



11 



